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We theoretically discuss one-photon and two-photon double-slit interferences for spontaneous and 
stimulated parametric down-conversions. We show that the two-photon sub-wavelength interference 
can exist in a general spontaneous parametric down-conversion (SPDC) for both type I and type II 
crystals. We propose an alternative way to observe sub-wavelength interference by a joint-intensity 
measurement which occurs for only type I crystal in a higher gain of SPDC. When a signal beam 
injects into the crystal, it may create two interference patterns by two stimulated down-converted 
beams, showing no sub-wavelength interference effect. 
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, I. INTRODUCTION 



The Young's doublc-slit interference experiment is one of the powerful ways to exhibit the nature of optical field, both 
^ ^ classical and nonclassical coherence effects. In the recent years, an interesting subject is devoted to the study of two- 
, photon double-slit interference in the process of spontaneous parametric down-conversion (SPDC). Q- Hi Since in this 
process a pair of converted beams created by a pump beam is in entanglement, the two-photon double-slit interference 
T-H may show some peculiar phenomena such as the quantum sub-wavelength lithography and the ghost interference. For 
• the former, both signal and idler beams are set together to pass through a double-slit, 0,0, 0-^3, ^"^^ fo'^ 
\ the latter, the double-slit is placed on the path for only one beam.0]-3,0],0|,E3 The original idea of the quantum 
<r lithography comes from the reduction of de Broglie wavelength for combining two massive particles. In optical system, 
the sub- wavelength interference occurs for a biphoton state. 16j In addition to the two-photon double-slit interference, 
Oh. the quantum lithography can be carried out in a Mach-Zehnder interferometer. Due to the fact that this effect 

' can overcome the Rayleigh diffraction limit, it may have prospective application in photo- lithography technology. 
^ \ In most of the theoretical analyses, the sub-wavelength interference is explained by a two-photon entangled state 
which can be acquired in the SPDC with very low gain. Nevertheless, it is the obstacle in practical application due 
to the lower power. Therefore, the exploration of these quantum effects in macroscopic regime makes sense. [l5 | . |18 | 
In this paper, we study one-photon and two-photon double-slit interferences in both spontaneous and stimulated 
parametric down-conversions. We focus on the case in which a double-slit is inserted on the paths for both signal 
and idler beams. We find that the sub- wavelength lithography can occur at very high gain of SPDC with substantial 
^ ' visibility. The discussion covers both type I and type II crystals which exhibit different behaviors in two-photon 
5^ interference. In the stimulated process, the amplified beam cannot perform quantum lithography but create rich 
interference patterns. The paper is organized as follows: in Sec. II we briefly review double-slit interferences for 
a coherent state and a two-photon state. In Sec. Ill we cite several formula as review of the optical parametric 
down-conversion process. We analyze two-photon double-slit interference in Sees. IV and V for the spontaneous and 
stimulated processes, respectively. The final section VI is the conclusion and discussion. 

II. DOUBLE-SLIT INTERFERENCES FOR A COHERENT STATE AND A TWO-PHOTON STATE 

We consider the scheme of Young's double-slit experiment as shown in Fig. 1. The double-slit function is defined 

by 

r(x) = rect(^-^) + rect(^±^), (1) 
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where d is the distance between the centers of two sUts and b is the width of each sht. In Fig. 1, both the double-sUt 
and the detection screen are placed at the two focus planes of a lens. By ignoring the thickness of the double-slit, the 
transverse envelope operators of the input field e{x,t) and the output field e'(x,t) of the double-slit are related as 

e'{x, t) = T{x)e{x, t) + [l- r(x)]e,,,(x, t), (2) 

where the vacuum field operator evac{x,t) is introduced for the sake of e'(x,t) satisfying the bosonic commutation 
relation. Since the vacuum field has no contribution to the normal-order correlation, it can be neglected in the 
calculations below. 

In the paraxial approximation, the field r{x,t) in the detection plane P2 is expressed by the Fourier transform of 
the lens 



r(a;,t) = y^^ J e'{x' ,t) exp[-ijx'x]dx' . (3) 

By substituting Eq. |(2J) into Eq. one obtains 

1 k f f ~ kx 



T{q) = / T{x)e-''i''dx = 4^ sinc(g6/2) cos{qd/2) (5) 

V 27r J V 27r 



where 



is the Fourier transform of the double-slit function T{x), and e(g, f2) is the Fourier transform of e(x,t) for both the 
spatial and temporal variables. 

First, we consider the input field to be a stationary, monochromatic plane wave in a coherent state 

{eix,t))^A, (6) 

where A is a constant. It has 

{e{q,n)) ^2TTAd{q)6{n). (7) 
In the detection plane, the first-order correlation is calculated as 

G(i)(xi,a;2,t) ^ (rt(xi,t)r(x2,t)) = M!f *(^)r(^). (8) 
Hence, the intensity distribution in the detection plane is written as 

G('\x,x,t) ^ {rHx,t)rix,t)) = ^f^^ = sinc^ ( ^ ) cos^ ( ^ ) (9) 

where /q = and A = 27r/fc. Note that f (a;) is a real function. Equation |(5J represents an interference fringe 

with the interval Xf /d in the range A//5, as shown in Fig. 2. 

Similarly, the second-order correlation function can be obtained as 



G'-^\xi,X2,t) = {rHxi,ty{x2,t)r{x2,t)r{xi,t)) 

- ^fH^)f\'^). (10) 
/ / / 

According to the theory of field coherence, the separability of spatial variables in the first-order and the second-order 
correlation functions verifies the perfect coherence of the field. Since the field operators at different positions are 
commutable, the second-order correlation of the field is in fact the spatial intensity correlation and it can be observed 
by a coincidence measurement as shown in Fig. la. The spatial patterns related to G^^\x,x,t) and G^'^\xi,X2,t) 
are called the one-photon and the two-photon interferences, respectively. According to Eq. (|10|) . in the coincidence 
measurement, if we scan one detector by fixing another, the same interference fringe as the one-photon interference 



3 



can be observed. Now, we introduce two other ways of observations for the two-photon double-slit interference. 
One is the spatial intensity-correlation measurement by scanning two detectors synchronously at a pair of symmetric 
positions, xi = —X2 = x. The other one is the two-photon intensity measurement by using a two-photon detector 
which generates an photo-electron by absorbing two photons. Applying these two observations to Eq. H1U|I . one has 

G(2) [x, X, t) = (x, -X, t) = II sinc4(^) cos'^(^). (11) 

In Fig. 2, we plot G^^^(x,x,t) (G*^^^ (x, — a;, i)) in comparison with G*^^-'(x,t) for the coherent beam. It shows that the 
two interference patterns are alike. 

The above discussion on the coherent state is analogous to the classical field. Then, we consider a two-photon state 
as input, which is the quantum state without classical analogue. A general two-photon state can be written as 

1^) = / dqsdq,C{qs,qi)al{qs)a\{qi)\{)), (12) 



where a\ and a\ are the creation operators for s and i photons which are assumed to be distinguishable, qs and qi 
are the transverse wave-vectors. When the input field is stationary, Eq. can be simplified as 



<^) = ^^ jT{j-q)eiq)dq. (13) 

By using Eq. p3(l , the first-order correlation functions for s-photon and i-photon in the detection plane are obtained 
as 

GW(a;i,X2) = dqdq,dq2C*{q,,q)C{q2,q)f*{'^ - q,)f{^ - q2), (14a) 

G['\x^,X2) = ^/ dqdq,dq2C*{q,q,)C{q,q2)f*{^ - q,)f{^ ~ q^), (14b) 

respectively. 

For two fields case, the second-order correlation function is defined by 

{xi,X2,t) = {r\{xi,t)rl{x2,t)rs{x2,t)r,{xi,t)^ , (15) 

which describes the coincidence probability of s-photon at position X2 and i-photon at position xi. In the case of 
Xi = X2 = X, it describes a two-photon intensity distribution. First, for the two-photon state H12I) . we calculate the 
two-photon wavepacket in the detection plane 

{0\rs{x2)n{xi)\'ilj) = — - / dqsdq,C{qs,qr)T{— qs)T{— q-i). (16) 

Then the second-order correlation is obtained as 

G^'\xi,X2) ^ mr,ix2Mx^m\^. (17) 

We discuss two extreme cases: two photons are independent and perfectly entangled in the transverse wave-vector. 
In the unentangled case, C{qs,qi) — C s{qs)Ci{qi), the first- and second-order correlations are written as 



G^''\xi,X2) = 



k 

k 
2^f 



, kxi 

T 



dqC:MT*C^-q) I dqC,n{q)T{ 



kx2 

T 



,kx2 



f 



, kx\ 



(to = s,i). 



dgG,(<z)T(— ^ - q) / dqCM)T{^ - ?) 



/ 



Gf^(a;i,:Ei)G«(x2, 



X2) 



(18a) 



(18b) 



Equations H18|l show two separabilities. The one is the separability of positions in both the first- and second-order 
correlation functions. The other is that the second-order correlation is factorized to two one-photon intensity dis- 
tribution. That is, the two-photon interference consists of two individual one-photon interferences, verifying the 
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Dirac's statement: "Each photon interferes only with itself. Interference between two different photons never occurs." 
However, these features mean the perfect first- and second-order coherencies for the two independent photons. 

In the opposite extreme, the perfect entanglement in wavevector, qs + qi = 0, is hold in state H12|l . For simplicity, 
we assume C{qs,qi) — > S{qs + qi), and the first- and second-order correlations are written as 



gW(xi,x.) = ar^(x,,x,) = ^|d,f*(^ + g)f(^ + ,) = -|=r[|(x2-xi)], 



(19a) 



X2) 



27r/ 



Jdqfi'-f+q)fi'-f-q) 



k ^2 r k 



27r/2 7 



r1-(xi+a;2)], (19b) 



where we use the integrals 



dqT*{^±q)T{^±q) = ^ dgdx;dx^r(a;;)r(4)e^(^±^)-i-(^±^)-^ (20) 



dx'idx'2T{x[)T(x'2)S{x'i ~ x'2)e^ 
= J dx{T{x[)T{x[)e'7^^'-^'^< = V2^f[j{x2 - xi)], 



and 



ldqf*i^±q)fi^Tq) = V2^f[^{ 



X2+Xi)]. (21) 



Note that T'^{x) = T{x). Equations (|19|l show that for the wavevector-entangled two-photon state the first- and 
second-order correlation functions in double-slit interference have position-correlation which results in decoherence. 
In the measurement, letting xi — X2 in Eq. H19|l . one obtains 

G(i) (x, x) = G« {x,x)^ , (22a) 

/v27r 



G^'\x,x) = ^T'[^i2x)]. (22b) 

Therefore, the one-photon double-slit interference disappears completely and the two-photon double-slit interference 
shows a sub- wavelength property since it has 

f^(|2.,c<si„c=|Jl|^l»^[J^|^|. (23, 

The fringe is the same as the ordinary double-slit interference with the half of the wavelength. The above analysis 
explains the complementarity of coherence and entanglement. 0,01,01 We emphasize that the discussion is also 
true for the case when s and i photons can be indistinguishable. 

III. THE BASIC FORMULA IN THE OPTICAL PARAMETRIC DOWN-CONVERSION 

In the optical parametric down-conversion, in which a plane-wave pum p field of frequency ujp activates a x'^' 
nonlinear crystal, the basic unitary transformation is described bv|l5j. poj-|22j| 

e„i{q,^) ^ U„i{q,n)a„i{q,n) + Vjn{q,^)al{-q,-^) {m^n = s,i), (24) 

where emil,^) and Umiq,^) are the output and input field operators, respectively, q is the transverse wavevector 
and is the frequency deviation from the carrier frequency. The transfer coefficients Umiq, ^) and Vm{q, ^) are given 

Us{q,n) ^ es{q,n)[coshr{q,n) + t^^^smhT{q,n)], (25) 

2T{q,n) 

V,iq,n)^e,iq,n)^,^smhTiq,n), (26) 
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u,{q, ii) = e,(g, r!)[coshr(-q, -n) + z^i^^ sinhr(-g, -n)], (27) 

V^iq, n) ^ 9,;(g, sinhr(-g, -f]), (28) 

where 



r(g, 17) = ^g^~A^q,n)/A, (30) 
A(q,l]) = [ks,{q,n) + h,{~q,--n) - kp%, (31) 

Ao = (fcs + ki- kp%. (32) 

g is the coupUng strength and Ic is the length of crystal. Aq is the coUinear phase mismatching of the central frequency 
components which correspond to the wave-numbers kj (j = s, i,p). For simplicity, we assume that two down-converted 
beams have the degenerate carrier frequency ujp/2. Hence, Eq. (|31|l can be reduced to an even function of both q and 

n 

Aiq,n)^ Ao + nynl-q^ql (33) 

where ilo and qo are defined as the frequency and spatial-frequency bandwidths, respectively. 

Equations (|24|I - H33|I describe the SPDC process of type II crystal but also can be applicable to type I crystal. 
For the former, two converted beams are orthogonally polarized, whereas for the latter, they are degenerate in both 
polarization and frequency but spatially separated. Therefore, Eq. (|24|l can describe type I crystal by omitting the 
subscripts. As a matter of fact, under the assumption of the carrier frequency degeneracy, Eqs. H25|) and (|26|) are the 
same as Eqs. and 



IV. DOUBLE-SLIT INTERFERENCE IN SPONTANEOUS PARAMETRIC DOWN-CONVERSION 

In Fig. 1, the two down-converted beams generated from the crystal illuminate a double-slit and then are detected 
in the focus plane of the lens. We designate 0^,(0:, i), emix,t), e'j^{x,t), rm(x,t) the slowly varying field operators for 
the input surface Pin and the output surface Pout of the crystal and the output plane of the double-slit Pi and the 
detection plane P2, respectively. Substituting Eq. (|24ll into Eq. Q), we may calculate the first- and the second-order 
correlations for the field in the detection plane P2. 

In this section, we consider the case of the spontaneous parametric down-conversion in which the input fields are 
in the vacuum state. The first-order correlations for two beams are obtained to be 

Gl^\xi,X2) = M„,{xi,X2) = {0\rl{xi,t)r^{x2,t)\0) {m = s,i) (34) 

= Wrll'^'"^^ \yr,Mm'T*i^~q)fi'^-q). 

Gm^ (x, x) illustrates the intensity distribution for beam m in the detection plane. Then we consider the second-order 
correlation function defined by Eq. H15|l . which now describes the intensity correlation between the signal beam at 
position X2 and the idler beam at position xi. In the case of xi = X2 = x, it describes a two-photon intensity 
distribution. For type II crystal, the two-photon intensity distribution can be measured experimentally by the scheme 
shown in Fig. lb in which the coincidence measurement of two orthogonally polarized photons is performed at the 
same position Xi = X2 = x. However, for type I crystal, the subscripts s and i should be omitted in Eq. (|15|) . If there 
is a two-photon detector available, one may observe the two-photon intensity distribution at position xi = X2 = x. 
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Alternatively, a realistic detective scheme for type I crystal is shown in Fig. la in which the intensity correlation is 
measured by two one-photon detectors at the different positions Xi and X2- 

For the vacuum state set in Eq. H15|l . by using Eqs. (0)) and H24|l . the second-order correlation is calculated as 

G^^\xi,X2) = M,{xuXi)M,{x2.X2) + \N,s{xi,X2f + 5,, \M{xuX2)\^ , (35) 
where 5is is 1 for type I and for type II crystal. Mm{xi,X2) is given by Eq. and 



N^n{xi,X2)^-^ JJ dqdn V,n{q,n)Unhq,~n)T{'^~q)T{'^ + q). {m^n = s,t) (36) 

Equation 1)35(1 shows that the second-order correlation is related to the first-order correlation. In Eq. ((35|l . the first 
term, which is separable in term of both polarization and position, describes the part of two-photon interference 
contributed by two individual single-photon double-slit processes. However, the second and the third terms describe 
two-photon interference effect related to photon entanglement. The difference between the two types of crystals is 
devoted by the third term, that is, the first-order correlation \M{xi,X2)\^- 

We note that, since Qi{q,V,)Qs{—q,—^) = Qs{q,^)'di{—q,—^) = exp[— iAq], it has Vi{q,il.)Us{—q,—il) ~ 
Vs{q^fl)Ui{—q,—fl), and hence Nis{xi, X2) = Nai{xi,X2)- This result is obvious because exchange of the indices 
i and s makes no difference physically. Moreover, due to Eq. (|33|l . A(q, i7) is the even function of q and under 
the assumption of the frequency degeneracy so that |t^i(q, f2)|^ = |t4(g, f2)|^. Therefore, two first-order correlation 
functions for the signal and the idler beams are equal in type II crystal, i.e. Mi{xi,X2) = Ms{xi^X2)- Nevertheless, 
we keep the subscripts in Eq. H35|) for a general description of type II crystal in case the two converted beams have 
different carrier frequencies. 

In order to obtain the analytical result for the integrals, we discuss two bandwidth limits of SPDC process: the 
broad and the narrow bandwidths. In the broadband limit, go 3> 27r/(i, Um{q,^) and Vmiq,^) are much flatter in 
comparison with T{q) and one can set Um{q, i^) ~ Um(0, ^) and Vm{q, f2) ~ Vm(0, il) in the integrals. By taking into 
account Eqs. (|20|l and (|21|) . Eqs. (|34|l and Ij^GI) can be rewritten as 

MUxi,X2) = ^Jdqf*{^-q)f{^-q) = ij„,f[j{x2-x,)], {m = s,t), (37a) 

and 



N.mn{xi,X2) ^ £,„inT[j{xi + X2)], {ui ^ n = s , i) , (38a) 

respectively, where we define rj^ = J |Kn(0, ^2)1^ dfl and C„„ = J Kn(0, ^)U„{0, -n)dn. 

In the broadband limit, which exhibits the maximum entanglement in transverse wavevector for two converted 
beam, we see again the position-correlation in the correlation functions. The first-order correlation Eq. H37a|) shows 
the same position-correlation as Eq. (|19a|l for the two-photon state with the maximum wavevector-entanglement. 
This makes the one-photon intensity distribution in the detection plane P2 homogeneous, Mm{x,x) — 77,„T(0), so 
that the one-photon double-slit interference disappears completely. 

In this limit, the second-order correlation (|35|l is obtained to be 



G(^)(xi,x2) = v,Vs {r'iO) + S^.T'[jix2 - xi)]} + \U' T'iji^^ + X2)]. (39) 

The first term in {} comes from two individual single-photon double-slit processes which are now homogeneous. The 
second term in {} and the last term manifest explicitly the position-corrlation. If we fix one detector at a position 
and scan another in the coincidence measurement, the interference fringe observed is the same as the single-photon 
one. To show the position correlation in the two-photon interference, we discuss two special observations, that is, 
xi = X2 = X and xi = —X2 = x. Equation H39|) is reduced to 

G^^\x, x) = 77.^.(1 + 5.s)f^{Q) + \ij f\j2x), (40) 

and 



G(2)(x,-x) = {r^.r^^ + \U')TH0) + 5.sri^VsT'i'^2x). (41) 
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The former exhibits a two-photon intensity distribution and the latter exhibits the intensity correlation at a pair 
of symmetric positions. Both Eqs. (|40|l and H41|l include a term T'^{j2x) which characterizes a sub-wavelength 
interference pattern by the factor of A/2 in comparison with the ordinary interference shown by Eq. Obviously, 
due to Eq. 1)40(1 , the sub- wavelength interference for the two-photon intensity distribution can occur in both type I and 
type II crystals. However, Eq. 1(41(1 shows that, for type I crystal, when a pair of single-photon detectors are placed 
at a pair of symmetric positions and moved synchronously in the opposite direction, the sub- wavelength interference 
can be also observed. But this effect never happens in type II crystal. 

According to Eqs. ((40(1 and 1(41(1 . the visibilities of fringes designated by G^^-'(a;,x) and G(2)(x,-a;) are calculated 
to be 



and 



- 1 + 2(1 + 

respectively, where 9 = rj^rjg/ IC^sl^. Note that V2 makes sense only for type I crystal, for which 9 = rf / |^|^. As the 
parameter 9 is increased from a very small quantity, Vi decreases monotonously from unity and V2 increases from zero 
up to 1/3. Since the parameter 9 is related to the coupling strength g of SPDC, we plot the visibilities as functions of 
g in Figs. 3, in which ViS for type II and type I crystals are indicated by the solid and dashed lines, respectively, and 
V2 for type I crystal is indicated by the dotted line. In a weak coupling of SPDC, which may generate approximately 
a two-photon entangled state, the visibilities Vi for both type I and II crystals reach perfectness. The sub- wavelength 
interference in the weak SPDC has been observed experimentally. l5l|,P|,'T^,'T^ The important fact is that, in Figs. 
3, the sub-wavelength interference can still exist, with a substantial visibility, even in very high gain of SPDC, in 
which the beams contain a large amount of photons. On the other hand, in the strong coupling of type I crystal, the 
sub- wavelength interference can occur in the joint-intensity measurement at a pair of symmetric positions. Equations 
((42(1 and ((43(1 show that, for type I crystal, two observations of sub-wavelength interferences compete for visibility. 
When Vi reaches the perfectness, V2 vanishes. However, in the very high gain of SPDC of type I, two visibilities are 
the same as 25%. 

In the opposite limit, we assume that the SPDC has a very narrow bandwidth go ^ 2Ti/d. Extremely, letting 
go ^ 0, the transfer coefficient Vm{q^ ^) {m = s,i) tends to the delta function 

Vr„{q,n)^ Vrn{0,n)S{q) (44) 

Equations 1(34(1 and ((36(1 are respectively written as 

M„.(xi,X2) = -^,,„f (^)f'^(^), (45) 

V 27r / / 

and 

iV_(xi,X2) = -L^_T(^)f(^). (46) 

V 27r / 7 

In this limit, the position-correlation disappears completely. The second-order correlation is then 

G(^)(xi,X2) = ^[(1 + <S,,)^,^, + \U']n'^)fH^). (47) 

Therefore, the one-photon intensity distribution M,„(a;,x) and the second-order correlation G^^-* (xi, 2:2) in the plane 
P2 are the same as the case for the coherent state. In the narrow bandwidth limit, since each converted beam is 
monochromatic, the one-photon double-slit interference occurs with the perfect visibility. On the other hand, two 
monochromatic converted beams have no more correlation in the transverse wavevector so that the position-correlation 
degrades completely in the two-photon interference. 

We plot the two-photon interference patterns by varying the bandwidth of SPDC process go in Figs. 4, in which 
Figs. 4a and 4b (4c and 4d) are for type I (type II) crystal. In Figs. 4a and 4c, a low gain of SPDC, g — (1/2) log 1.5 
(with the amplification rate 1.5) is taken in two-photon intensity measurement so that the sub- wavelength interferences 
with a better visibility are achieved when the normalized bandwidth qgh /{2Tr) is increased. Two sets of patterns are 
very similar with exception of tiny higher intensity for type I crystal. Figures 4b and 4d show the interference patterns 
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for the joint-intensity measurement of two one-photon detectors at a pair symmetric positions. The sub-wavelength 
interferences can be observed only for type I crystal when the gain of SPDC is higher, for instance, g = (1/2) log 10 
(with the amplification rate 10) is taken in the figures. Though the visibilities are lower, the intensities of the patterns 
are getting much higher. This also happens in the case of two-photon intensity measurement. The three plots of Figs. 
4a-4c show that the bandwidth of SPDC governs two-photon sub- wavelength interference. For a very small bandwidth 
of SPDC, the two converted beams are de-entangled in transverse wavevector so that the nonclassical sub- wavelength 
interference disappears. 



V. DOUBLE-SLIT INTERFERENCE IN STIMULATED PARAMETRIC DOWN-CONVERSION 



In the stimulated optical parametric process, a signal beam is injected into the nonlinear crystal and is then 
amplified. The nonlinear crystal becomes an optical parametric amplifier (OPA). We assume inputting a stable 
plane-wave beam in a coherent state 

{Zs{q, n)) = 2TrAS{q ~ Q)S{n), (48) 

where Q designates the transverse wavevector of the input beam deviated from the normal incidence. For type II 
crystal, we set the input beam as the signal which can be identified by the polarization, while the idler beam is in 
the vacuum state. For type I crystal, the subscript s in the input beam (|48(l is omitted. Considering the input beam 
describing by Eq. (|48|) instead of the vacuum state, we calculate the first-order correlation in the plane P2 



where 



for type II crystal, and 



Gl^\xi,X2)^W:;^ixi,Q)Wmix2,Q) + M,^{xi,X2), (m = s,z) (49) 



W,ix.,Q) = A^/k/JUs{Q,0)T{kx/f-Q) (50a) 
W,{x,Q) = A^k/Jv^i-Q.O)f{kx/f + Q) (50b) 



W{x, Q) = A./k/J[U{Q, 0)T{kx/f - g) + V{^Q, 0)T{kx/f + Q)] (51) 



for type I crystal. Mm{xi,X2) has been defined by Eq. In Eq. the first and the second terms show 

contributions coming from the stimulated and the spontaneous processes, respectively, and they are independent 
in the onc-photon interference pattern. Obviously, the stimulated part shows the first-order coherence due to the 
separability of the spatial variables. Nevertheless, for the input beam is enough strong, the spontaneous process 
can be neglected. Gm\x,x) ~ |Wm(a;, Q)P describes an amplified double-slit interference pattern in comparison 
with the case without crystal. For type II crystal, the input signal beam creates two interference patterns: one for 
the signal beam with the amplification rate \Us{QjO)\'^ and the other for the idler beam with the amplification rate 
\Vs{—Q, 0)p.|23] According to Eqs. H50|l . these two patterns are the same as the ordinary fringe (see Eq. @) and can 
be identified by the polarization and separated in space when the input beam is well tilted in incidence. However, for 
type I crystal, the one-photon interference pattern is written as 

\W{x,Q)\' = !l±{\U{Q,0)\'T\j-Q) + \V{~QM'T'{j + Q) 

+ [U{Q, Q)V*{-Q,0)f{j - Q)f{j + Q)+ c.c.]}. (52) 

The first and second terms illustrate the two interference patterns created by two stimulated beams, i.e. the signal 
beam with the transverse wavevector Q and the idler beam with the transverse wavevector —Q. The third term 
shows an additional coherent superposition of the two interferences of the converted beams. In result, the interference 
pattern can be different from the ordinary one due to the " interference term" in the square brackets of Eq. (|52(l . Only 
for the normal incidence, (5 = 0, the interference pattern is the same as the ordinary one. Figures 5a and 5b show 
the density patterns of the one-photon interference by varying the transverse wavevector of the input beam for type 
I and type II crystals, respectively. In Fig. 5a, the fringes alternate between onset and offset when the transverse 
wavevector Q of the input field is increased until the two fringes are well apart. The offset is due to the destructive 
interference of two indistinguishable stimulated beams when they are folded. However, for type II crystal, we put two 
interference patterns for the signal and idler beams together in Fig. 5b, in which the right part is for the signal field. 
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and is stronger than the left part (idler). This corresponds to the measurement when the detection is insensitive to 
the polarizations. Different from type I crystal, when the transverse wavevector Q of the input field is increased, the 
bright and dark spots of the patterns exchange alternatively until the two fringes are apart. This feature comes from 
the incoherent superposition of the two patterns for the signal and idler beams. 

We go through a long derivation, using the unitary transformation (|24|l and the bosonic commutation relation, and 
obtain the second-order correlation function 

G(2) {xi,X2) ^[\W,{xi,Q)\^ + Nh{xi,xi)][\Ws{x2,Q)\^ + Ms{x2,X2)] (53) 

+ [W^ {xi,Q)Ws {X2 ,Q)Nl{xi,X2)+C.C.] + \N,s{xi,X2)t 

+ 6,,{[W*{xi,Q)W{x2M)M{xi,X2) + C.C.] + \M{xi,X2)\'}. 

Again, Eq. (|53|l can describe two types of crystals. For type I crystal, the subscripts i and s are omitted and 
Sis — 1, whereas for type II crystal Sis = 0. Similar to Eq. H35() . the first term is the product of the two one-photon 
interferences for the converted beams. The second and the fourth terms exhibit the coupling of the interferences 
between the stimulated and spontaneous processes. This result is similar to the discussion of the image amplification 
in the optical parametric amplification in which W{x) describes an amplified image. p3| We note that the interferences 
coming from the stimulated process are uncorrected in positions. This reflects that the stimulated process does not 
includes the entanglement in transverse wavevectors for the two converted beams. Therefore, the quantum lithography 
cannot be achieved in the stimulated parametric amplification by input a coherent beam. In Figs. 6, we plot the 
density patterns of the two-photon interferences by varying the input direction of the signal field. Figures 6a and 
6b show patterns of G'^^ (x, x) and G'^^ (x, —x) for type I crystal, respectively. These two patterns are similar to 
the one-photon interference case with the exception of that the sub-wavelength fringe coming from the spontaneous 
process appears in the central part. In Fig. 6a, the fringe of the signal part (the right side) is stronger than the idler 
part (the left side), whereas in Fig. 6b, the pattern is mirror-symmetry since exchange of two detectors at a pair of 
symmetric positions makes no difference. However, Figs. 6c and 6d show patterns of G^'^\x,x) and G(^^(a;, —x) for 
type II crystal, respectively. Correspondingly, Fig. 6c is similar to the one-photon interference case with the exception 
of that of the spontaneous contribution. In Fig. 6d, we define x as the transverse position for the detector sensitive 
to the polarization of the signal beam, so that it is comparable to the one-photon case for the signal beam (the right 
part in Fig. 5b). Furthermore, in Fig. 6d, the spontaneous process contributes a homogeneous background instead 
of a sub-wavelength fringe in Fig. 6c. If two detectors in the joint-intensity measurement are polarization-insensitive, 
one observes a symmetric patterns which is the mirror-symmetry of Fig. 6d. 

VI. CONCLUSIONS 

In summary, we formulate the first- and second-order correlation functions in the Young's double-slit interference for 
both spontaneous and stimulated parametric down-conversions. The results reveal the relations between the first- and 
second-order correlations, and hence it can explain the complementarity of coherence and entanglement. We show 
that the nonclassical sub-wavelength two-photon interference can occur macroscopically in a general spontaneous 
parametric process. For a high gain of SPDC, in which the converted beams contain a huge number of photons, the 
sub-wavelength interference pattern is intensive with a substantial visibility. This makes the quantum lithography 
technology in practicability. Moreover, we find an alternative way to observe the sub-wavelength interference for 
only type I crystal, in which a joint-intensity measurement is performed by a pair of one-photon detectors placed 
at the symmetric positions. The advantage of this method is that the quantum lithography in type I of SPDC can 
be performed by two one-photon detectors instead of a two-photon detector which could be unavailable. Since, this 
effect occurs only in a higher gain of SPDC, it reflects macroscopically quantum nature for the entangled beams 
containing a huge number of photons. The two ways of observations compete for the visibility, so that at the low gain 
of SPDC of type I the interference of the first observation reaches perfectness while the second observation disappears. 
In the stimulated process, the one-photon and two-photon interference patterns generated by a pair of stimulated 
down-converted beams are alike. For type I crystal, since two converted beams are polarization-indistinguishable, the 
coherent superposition of two converted beams causes a secondary interference which may fade the fringe when two 
beams are not well apart. However, this effect does not exist for type II crystal because of the distinguishability in 
polarization for the two stimulated beams. 
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Captions of Figures : 

Fig. 1. Schemes of Young's double-slit interference with a convex lens: (a) a one- photon (two- photon) detector 
measures one-photon (two-photon) intensity distribution; two one-photon detectors measure joint-intensity distribu- 
tion at a pair of symmetric positions, (b) for type II crystal, two-photon intensity distribution is measured by a 
polarizing beamsplitter (PBS) and two one-photon detectors. 

Fig. 2. One-photon (solid line) and two-photon (dashed line) doublc-slit interference patterns for a coherent beam. 

Fig. 3. Visibilities versus the gain of SPDC for different coUinear phase-mismatching (a) Aq = —5.85; (b) Aq — 0; 
and (c) Aq = 5.85. Solid and Dashed lines designate Vi for type II and type I crystals, respectively; Dotted line 
designates V2 for type I crystal. 

Fig. 4. Two-photon interference patterns versus the normalized bandwidth of SPDC qob/{2Tr): (a) G'-^\X,X) for 
type I crystal; (b) G(2)(A,-A) for type I crystal; (c) G^'^^X.X) for type II crystal; and (d) G^'^^X,-X) for type 
II crystal, where X — xkb/{2Trf) is the normalized transverse position in the detection plane. The gains are set as 
g = (1/2) log 1.5 in (a) and (c), and g = (1/2) log 10 in (b) and (d). The other parameters in Figs. 4-6 are taken as 
the phase matching Aq = 0, the ratio b/d — 0.2, and = 1 for an arbitrary unit of the correlation functions. 

Fig. 5. Stimulated one-photon interference patterns by varying the normalized transverse wavevector of the input 
field Qb/{2'K) for (a) type I crystal and (b) type II crystal. The normalized bandwidth and the input intensity are 
taken as qob/{2TT) — 2 and ^^A"^ — 1, respectively. 

Fig. 6. Stimulated two-photon interference patterns by varying the normalized transverse wavevector of the input 
field Q6/(27r): (a) G^'^\x,x) for type I crystal; (b) G^'^\x,—x) for type I crystal; (c) G^'^'>{x,x) for type II crystal; 
and (d) G'^^-'(x, —x) for type II crystal. The normalized bandwidth and the input intensity are the same as in Fig. 5. 
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